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Abstract 



The zero-temperature system of two dilute overlapping Bose-Einstein con- 
densates is unstable against long wavelength excitations if the interaction 
strength between the distinguishable bosons exceeds the geometric mean of 
the like-boson interaction strengths. If the condensates attract each other, the 
instability is similar to the instability of the negative scattering length con- 
densates. If the condensates repel, they separate spatially into condensates 
of equal pressure. We estimate the boundary size, surface tension and energy 
of the phase separated condensate system and we discuss the implications for 

double condensates in atomic traps. 
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As dilute gases, the atomic trap Bose-Einstein condensates occupy a unique position 
among the superfluid systems. One intriguing consequence of their dilute gas nature is 
the prospect of studying condensate mixtures, which are the first experimentally realizable 
bosonic superfluid mixtures 0-0 (the only other mixture of superfluids is the fermion- 
boson 3 He- 4 He system). Understandably, this prospect has attracted interest J7J- [ID|, and 
recently the first observation of overlapping condensates was reported f|. 

In this paper, we study if and when zero temperature dilute condensates overlap. We find 
that the homogeneous overlapping condensate system is unstable against long wavelength 
excitations if the strength of the interaction between the distinguishable bosons exceeds 
the geometric mean of the like-boson interaction strengths. In that case, two repelling 
condensates spatially separate into single condensates of equal pressure. The condensates 
still partially overlap in the boundary region that separates them. We estimate the size of the 
boundary region, as well as the corresponding surface tension. We show how the homogenous 
treatment may be generalized to describe phase-separated large double condensate systems 
in a trap. 

The instability of the overlapping condensates manifests itself in the energy dispersion 
of the elementary excitations |10|| , as we show below. The wavefunctions, 0i and 02, of two 



interacting condensates satisfy coupled Gross-Pitaevski equations: 
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where fij (j = 1, 2) represents the chemical potential of the j-bosons. The interaction 
strength values, Xj and A, are determined by the scattering lengths for binary collisions of 
distinguishable bosons: Xj = 4irh 2 aj/mj and A = 2irh 2 a/m re d, where m~} d = m^ 1 + m^ 1 



The excitations of the static homogeneous condensates, <j)j(r,t) = <pf\ are described by 
fluctuations of the fields, <pj(r,t) = <f)f^ + 5<pj, which evolve according to the Gross-Pitaevski 
equations (JJ), linearized in 5(f) and 5<p* [|TTJ. Decomposing the field fluctuations into Fourier 



components, = J2k c j,k exp(zk-r), we obtain the equations of motion for the c-amplitudes, 
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where n, 



6^1 2 and where we have used that 0° 



0. A second equation for ihc 2 is 



obtained by exchanging the 1 and 2 subscripts. Alternatively, we can introduce the phase 
and density of the condensate field, <p = v /pexp(i6 ) ), the fluctuations of which, p = n + 5p 
and 9 = 0(°) + 59, account for the field fluctuations, 5<p = (o) [Sp/2n + i59). The last 



term of Eq.(El) then represents a density fluctuation, 5p k 



'c k + 



whereas the 



difference, <5n k = 0(°)*c k — <p^c*_ k /2i, represents a phase fluctuation. Multiplying Eq. 



by 



and adding and subtracting the resulting equalities with the complex conjugate 



equation, the 511 and dp -equations of motion follow: 
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Thus, the phase fluctuations of one condensate couple to the density fluctuations of the 
other. We cancel out the dependence on the phase fluctuations by taking the derivative of 
the first equation in Eqs. (|]) with respect to time and by substituting 5Il k from the second 
equation. With 5pk(t) = <5p k cos(f2 k t), we find the normal mode equations for the coupled 
density fluctuations, 
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where fwjj^ = \J \h 2 k 2 /2mj) 2 + (h 2 k 2 jm^)n^\j (we assume Aj > 0) denotes the usual single 
condensate Bogoliubov dispersion. Requiring Eq.(f§) to have non-trivial solutions gives the 
dispersions of the double-condensate excitations: 



where Cj is the sound velocity of the j-condensate, Cj = JrijXj/mj. The ±-sign in the sub- 
script, Q±, corresponds to the choice of + or — in Eq.(||). The implications for the physics of 
the overlapping condensate systems are profound: the near-equilibrium dynamics and ther- 
modynamics of the two-condensate system are essentially the dynamics and thermodynamics 
of two sets of nearly non-interacting quasi-particles with energies fl±k- 

We are concerned with the stability of the overlapping condensate system to which 
purpose we consider the long wavelength limit, k — > 0, of the fi-dispersion relations. With 
the long wavelength Bogoliubov energies, Uj^ ~ Cjq, we find that the double condensate 
dispersions of Eq.(Q) are also phonon-like, f2 ± k ps c±k (k — > 0), with 'sound velocities' c + 
and c_ determined by 



c ± = 2 • (6) 

If A 2 > A1A2, the inter-condensate interaction repels the energy levels so strongly that c 2 
becomes negative. Studying the normal mode frequencies of small deviations is a standard 
test of stability and c 2 < indicates that the overlapping condensate system is unstable 
against long wavelength excitations. 

The thermodynamic properties of the double condensate are determined by minimizing 
the free energy F. If the spatial variations of the condensates are slow and the kinetic energy 
contributions may be neglected, the zero-temperature free energy is the integral over the 
free energy density F(r), 

F(t) = Y n i( r ) + y n 2( r ) + ^i(r)n 2 (r) 

-f/i(r)m(r) - /i 2 (r)n 2 (r) , (7) 

where the effective chemical potentials fij(r) include the external potentials Vj(r) experienced 
by the j - bosons, /Xj(r) = Hj—Vj{r). Minimizing F with respect to the densities, 6F/5rij(r) = 
0, gives the Thomas-Fermi equations, 

//i(r) = Aini(r) + An 2 (r) , 

A* 2 (r) = X 2 n 2 (r) + An x (r) . (8) 



When Vj(r) = 0, the Thomas- Fermi condensate densities of Eq.(j8|) are homogeneous. How- 
ever, equating first-order derivatives to zero, only gives a minimum provided the second-order 
derivatives satisfy (d 2 F/dn 2 ) > and (d 2 F / dn 2 )(d 2 F / dn 2 ) - (d 2 F / dn^) 2 > 0. The lat- 
ter condition implies that the Thomas-Fermi equations (^j) only gives a minimum provided 
the stability criterion, A 2 < A]A 2 , is satisfied. 

To see that 'strongly' repulsive condensates, A > \/A]A 2 , lower their free energy by 
distributing the condensates inhomogeneously, we write the free energy density of Eq. (0) in 
the absence of external potentials, Vj(r) = 0, as 



Starting from the homogeneous overlapping condensate system, redistributing bosons 1 and 



2 spatially while keeping [^i(r) + n 2 (r) J A2/A1] constant over space, can lower the energy 
by decreasing the overlap integral / d 3 r ni(r)n 2 (r). The lowest value is reached by spatially 
separating the two condensates. In that case, although the hamiltonian is translationally 
invariant (vj(r) = 0), the double condensate is not: the double condensate system sponta- 
neously breaks translational symmetry. Note that the phase separation is consistent with 
the above discussed dynamical instability - the separated condensate system has no region 
that can be described locally as two homogeneous condensates. In contrast, strongly at- 
tractive condensates (A < 0, A 2 > AiA 2 ) decrease the free energy by increasing their mutual 
overlap and overlapping condensates are unstable. This is similar to the behavior of single 
condensates of negative scattering length, and we expect the strongly attractively double 
condensates to be similarly unstable. The same analogy suggests that a confining potential 
might give stable or metastable strongly attractive double condensates. 

The free energy of the phase separated condensate system in a macroscopic volume V, 
is the sum of the single condensate free energies, condensate 1 confined to a volume V%, and 
condensate 2 confined to V — V\. Minimizing the total free energy with respect to V\ gives 





(9) 
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the equilibrium condition of equal pressures exerted by both condensates. With the pressure 
Pj = \jn? s /2, where n^ s denotes the density within the separated condensates, we are lead 
to the equivalent condition for the condensate densities ni iS = n^sv/^/Ai. 

We note that the assumption of slowly varying condensate wavef unctions, necessary to 
justify neglecting the kinetic energy in the free energy of Eq.(|7]), is violated at the boundary 
of the two condensates. In fact, an infinitely sharp boundary gives an infinite kinetic energy 
contribution. The effect of the kinetic energy is then to give a boundary region of finite 
size b, in which the wavefunctions smoothly tend to zero as the condensates cross the 
boundary. We assume that the condensates are so large that the boundary region which 
separates them can be approximated locally as a planar region with densities that vary 
spatially as functions of the coordinate z with the z-axis perpendicular to the boundary 
surface of area A. The kinetic energy contribution, E^ in (b\ is then approximately equal to 
Ekinip) « (Ah 2 n ltS /2mib)[l + y / Ai/A 2 (mi/m 2 )], where we used that (n 2)S /ni )S ) = y^Ai/%. 
The overlap of the condensates in the boundary region increases the interaction energy by 
an amount E int (b), which we estimate by modeling the condensate densities in the boundary 
region, z G (0,6), crudely as n%(z) ~ n\ >s {b — z)/b and n 2 (z) ~ n^ s z/b. With Eq.(Q), we 
find that Ei nt (b) ~ A(X — VAiA 2 )ni jS n 2iS 6/6. To estimate the actual boundary size b, we 
minimize the boundary energy, Ej,(b) = Eki n {b) + Ei nt (b) with respect to b and find 

b = 2hV3 

where l± is the coherence length of condensate 1, l\ = h/ J4m\n\ iS \\. The boundary con- 
tribution to the energy is a surface energy E b (b) = a A, where the surface tension a is 
proportional to the coherence length /i, to the pressure P±, Pi = n\ s \x/2 y and to a di- 
mensionless constant Si which depends solely on the mass and interaction strength ratios, 
Si = 401 + (m 1 /m 2 )^Ai/A 2 ][A/ v / A^- lJ/v/3, o = hP^. 

In the absence of external potentials, a 'droplet' of condensate 1 immersed in a much 
larger condensate 2 of density n 2 minimizes Ej, by taking on the shape of a sphere of radius 
R, V\ = (An /3)R 3 . We can now imagine creating the double condensate system starting 
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from a single condensate of type 2 and replacing condensate 2 bosons in the droplet volume 
V\ by condensate 1 bosons. The energy AE required in the replacement is equal to AE = 
[Ai jS nf/2 — A2?i| s /2]Vi + Eb. Minimizing the 'replacement energy' AE with respect to V\ 
and realizing that Eb oc V^ 3 , we find 

Aiwf ia = \ 2 n\ s 2 E b 

2 2 3 Vl ' [ ' 

The previous result, \in\ s /2 = \2n\gj2, obtained by ignoring the boundary energy, is ac- 
curate provided the size of the droplet exceeds R s = 2a /Pi = 2l{Ei. The energy per droplet 
particle, AE/N ly with Eq.flTTD is equal to AE/Ni = 5E b /Vi = 5a/[Rn ltS ], a function that 
decreases monotonically as Ni increases. Consequently, splitting up the droplet into smaller 
droplets further increases the free energy and it is energetically favorable for condensate 1 
to gather in a single region of space (i.e. real space condensation). 

To describe separated double condensates in traps, we subtract the overlap term, Ani?i2, 
and include the boundary surface energy in the expression of the free energy (Eq.(^)). The 
validity of this description rests on two conditions: 1. the local coherence length within 
each condensate is much less than the length scale on which the condensates vary spatially 
and 2. the change of the potential energy across the inter-condensate boundary, \fj\b, where 
fj = —Vvj represents the external force experienced by bosons j near the boundary, is 
much less than the local chemical potential \ij\b << Xjrij S (j = 1, 2). If these conditions are 
satisfied, the physics of the phase separation is similar to the above v j(r) = 0-case, and we can 
answer interesting questions regarding trapped phase separated condensates. For instance: if 
we add a droplet of condensate 1 to a trapped condensate 2, does it 'sink' to the middle of the 
trap, or does it remain 'floating' on the surface of condensate 2? For the sake of simplicity, 
we assume that the size of the droplet is large enough to neglect the boundary surface energy 
and small enough to neglect the spatial variation of the density inside the droplet. Then, 
the previously defined 'replacement energy' AE depends on the center of mass position R 
of the droplet through the external potentials, fi(r) = v (r), and ^(r) = av(r). Since the 
pressures inside and outside the droplet are equal, Ain^ 2 /2 — X 2 n\ s /2 w 0, we find that 

7 



AE(R) = f [m(r)vi(r) - n 2 (r)t> 2 (r)] d 3 r 



^(^[l-aJVAa] • (12) 



Thus, even though both bosons experience a trapping potential, if aJ Ai/A 2 > 1, the force 



on the droplet, — A^fl — ayAi/A 2 ]Vv, is directed outwards : the droplet 'floats'. 

The thinner the layer of the floating condensate, the more important are the effects 
of the boundary surface energy. Indeed, with a few condensate 1 particles, the energy 
will be minimized by covering only part of the surface of condensate 2 (thereby reducing 
the surface boundary energy). Here, we only consider the case where enough bosons have 
been added for condensate 1 to 'wrap' around condenate 2, and we can ignore boundary 
surface energy effects. Within the single condensate regions, the condensates are described 
in the Thomas-Fermi approximation of Eq.(|) (putting A = 0), ni(r) = [/ii — t>(r)]/Ai 
and n 2 (r) = a\\i 2 — v (r)]/A 2 , where /ii and a\i 2 are the chemical potentials. In that case, 
the boundary surface is the equipotential surface v (R) = /if, where the pressures of both 
condensates, Ajn|(r)/2, are equal. This leads to 

\fii - /x 2 ] 



(13) 



ayAi/% - 1 

In figure 1, we show a typical density profile for two separated condensates in a spherically 
symmetric trap. With Xjnj(Rb) = fij — Vj(Rt), we find for the case shown in Fig.(l) that 



A*i — A*2 — ^i n i(Rb) — A 2 ra 2 (i? fe )/Q; > 0. Using [ni(R b ) / n 2 (Rb)] = J A 2 /Ai we then find that 



ayAi/A 2 > 1, so that condensate 1 should indeed float on top of condensate 2. Of course, 
the experimentally relevant quantities are the number of boson particles, Ni and N 2 , rather 
than Hi and /i 2 . The chemical potentials can be determined by inverting Ni = J Vi d 3 r ni(r) 
and N 2 = fy 2 d 3 r n 2 (r), where the boundary between V\ and V 2 is defined by Eq. fli"3]) , to 
give Nji/j,!,^)- 

In summary, we have shown that the homogeneous overlapping double condensate system 
with strong inter-condensate interactions (A 2 > AiA 2 ) is unstable. In that case, attractive 
condensates (A < 0) collapse and repelling condensates (A > 0) separate spatially. In 
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between the separated condensates is a region of partial overlap. We have estimated the size 
of this region, as well as the resulting surface tension. Finally, we briefly discussed phase 
separation of large condensates in atomic traps. 

The phase separation suggests many experimental applications. For instance, one can 
use a condensate to spatially confine a droplet of a different condensate in the middle of 
a trap. Using light that is resonant with the droplet atoms, the droplet can be displaced 
and its subsequent motion inside the confining condensate can be observed. If the motion 
is undamped, we have a direct observation of superfluidity. To further motivate such ex- 
periments, we mention that damping generally could occur at velocities less than the sound 
velocity of the confining condensate because of the creation of superfluid vortices etc... Fur- 
thermore, it is interesting to note that schemes have been proposed to continuously alter the 



interaction strengths, either by using light ||13|| , or by varying the bias field of the magnetic 
traps |L4] , so that it might be possible in the future to observe phase separation in real time 
as the interaction strengths are altered. 
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Figure Captions 

Fig.l : Plot of a typical phase separated double condensate in a spherically symmetric trap. 
In reality, the boundary of condensates 1 and 2 (at R — R^ where v(Rf,) = fib) is not 
infinitely sharp and the condensates overlap over a region of size b that is estimated in the 
text. 
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